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Abstract 

It is recalled why, as already stated in a previous paper, there seems 
to be an inconsistency in identifying the non perturbative Zf OM as 
the renormalisation of the axial current, or equivalently, in setting as 
normalisation condition that the renormalised vertex = 1 at p 2 = g 2 
at some renormalisation scale p, where p is the momentum in the legs. 
Indeed, unlike the vector case, the Ward-Takahashi (WT) identity for 
the axial current is shown to imply both the renormalisation scale 
independence of Za and a p 2 dependence of Z^ OAI . This p 2 depen¬ 
dence is simply related to certain invariants in the pseudoscalar vertex 
and can persist in the chiral limit due to the spontaneous breaking of 
chiral symmetry (pion pole). It is seen clearly in the p 2 dependence 
of some lattice calculations of Z^ OM /Zy OM near the chiral limit. 
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1 Introduction 


The non perturbative MOM renormalisation scheme for lattice as introduced 
in Martinelli et al. Nucl.Phys. B445 (1995) 81-108, [lj, and inspired by the 
corresponding continuum scheme of Georgi and Politzer [9], has represented 
an imprtant progress in lattice QCD. It is very intuitive and easy to handle by 
pedestrians because it relies on the vertices and propagator functions which 
are familiar in pertrubative QCD. This is to be compared with later schemes 
like the one of Alpha [2], which is rigourous, but requires much more effort 
to understand and many technicalities. 

Nonetheless, it has features coming from its direct hadronic meaning that 
precisely complicate the matching with the usual perturbative schemes which 
is the final goal. It has a non logarithmic, power dependence on yii 2 generated 
by OPE power corrections (4| and particle poles (although the latter are 
lying outside the Euclidean range). It may have critical chiral behaviour 
unlike the Alpha prescription, as has been underlined some time ago for the 
pseudoscalar Z P [3], leading to the recipe of ’’extracting the pion pole”. 

Moreover, it seems to have often escaped the attention that there is also in 
principle an inconsistency in introducing Z^ IOAI as the renormalisa¬ 
tion of the axial current, as stated in our paper PhysRevD.81.094504 [5j. 
This manifests itself in a /i 2 dependence of Z^ rOM /Zy OM which persists in 
the chiral limit, even with an explicitly chiral invariant action, due to the 
spontaneous breaking of the symmetry, see our earlier paper [6j, section 8 . 1 .. 

Since this statement was presented in a rather long paper devoted to 
several topics, we think useful to recall the arguments in a clearer and more 
explicit manner, and in the same course to correct some sloppy notations of 
the paper. 

It must also be said that having reread the basic paper [I], we have 
rediscovered that there was a discussion in it having connection with the 
present one, although the conclusion seems different : we extract a finite 
effect in the chiral limit, which does not appear in their approach. The 
reason will appear after having presented our own discussion, in a separate 
section 0 It requires an examination of the interplay of the chiral and q —> 0 
limits. 


2 Practical meaning of the problem 

One must warn from the beginning about the practical meaning of the prob¬ 
lem . The problem is found to disappear at large /x 2 , so that one may claim 
that it is not real since anyway, the non perturbative MOM scheme is meant 
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precisely to be applied at such large p 2 . More precisely, for Zy,A one seems 
to be free to choose any p 2 , therefore, it would suffice to work at such a large 
/i 2 with Z^ OM . In fact the effect seems anyway to be small already about 
p = 2 GeV or beyond, see fig. 5 in the quoted paper |6j, which is now the 
commonly adopted value for non perturbative renormalisation. Therefore, 
the effect is perhaps not worrying practically, in contrast to the pion pole in 

yMOM 

However it is worth in general being clear on theoretical principles. But 
also working at large p 2 presents well-known practical problems. 

1) one may always wonder how large p 2 must be. 

2) measuring Green functions at large momenta requires the extraction 
of large artefacts, and this was the initial motivation for non perturbative 
renormalisation : to avoid working at too large momenta. 

On the other hand, if p 2 is not sufficiently large, one has to “extract” 
a physical /i 2 dependence, while artefacts may still be non negligible (one 
cannot exclude non canonical artefacts at small /i 2 ). 

For all these reasons, it is useful to know about the possible causes of 
momentum dependence, either “physical” as the present one (physical with 
many quotation marks), or artefactic. 

One must add that at present little effort has been devoted to determine 
the actual magnitude of the effect (see the end of the text). 

3 Definitions and generalities 

Let us first fix the notations that we will use. We will use all along the 
Euclidean metrics. The continuum quark propagator is a 12 x 12 matrix 
S(pf,) for 3-color and 4-spinor indices. One can take into account Lorentz (in 
fact 0(4)) invariance and discrete symmetries, as well as color neutrality of 
the vacuum by expanding the inverse propagator according to : 

S- 1 (p) = S a _bZ t (p 2 )(i^ + M(p 2 )) (1) 

where a, b are the color indices. “Z^(p 2 )” is a standard lattice notation, 
alluding to the role it plays as a renormalisation constant for the quark field 
in the standard Georgi-Politzer MOM renormalisation, where Z| /0iU (/i 2 ) = 
Z^/r 2 )- 1 . But let us stress that here it is not by itself a renormalisation 
constant. On the other hand, M(p 2 ) is the mass function, which is one 
possible concept of mass, introduced by Georgi and Politzer. M(p 2 ) is UV 
finite, since it is the ratio of two quantities renormalised by the same factor 
Z 2 . It is identical with the MOM renormalised mass at scale p 2 = p 2 , see 
also below. 
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Let us consider a colorless local two quark operator qOq. The correspond¬ 
ing three point Green function G is defined by 


G(p, q ) 


d 4 xd 4 y e ip ' y+iq ' x < q(y)q(x)Gq(x)q( 0) > 


( 2 ) 


It is a 4 x 4 matrix in Dirac space. The associated vertex function is then 
defined by amputation of quark propagators on both sides : 

T (p,q) = S-\p)G(p,q)S-\p + q) (3) 


In the note, we will often restrict ourselves to the case where the operator 
carries a vanishing momentum transfer q^ = 0. We will then omit to write 
q^ = 0 and we will moreover understand T(p) without R index as the bare 
vertex function (computed on the lattice). 

Now, Lorentz covariance and discrete symmetries allow to write for the 
axial vertex r Al Jp, q) at q — 0: 

Ta/,(p) = ^a,b[9 { A (p 2 h^5 + ig { A\p 2 )P ^75 + 

9 { a O 2 )Prihs + ig { A (p 2 ) [ 7 #*, i>\ 75] (4) 


which should be obeyed approximately on the lattice, as we checked. 
On the other hand, we need the pseudoscalar vertex at q ^ 0 : 


r 5 (p,q) = S a ,b 


q )75 + ig {2 \Pi <?)75(Pm7/J + 


wt\p, q) ihiyft + gt\pi q)iAi^^\ 


(5) 


where the quark momenta are p,p + q and the g$\p, q)’s are invariant func¬ 
tions of the momenta alone. For brevity, the first two invariants are denoted 
by the same symbol at q — 0, i.e. : 


gi 1,2 \p 2 ) = g { 5 L2 \p,q = o) 

More explicitly, the dependence of (p, q) in p, q is : 

gfip,g) = ^(pW.p-s) 


( 6 ) 


(7) 


In the following, the color factors 5 a ,b will be skipped. 
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3.1 Renormalisation in general 

Without requiring any specific renormalisation scheme, we have to refer to 
the renormalisation, because the Ward-Takahashi(W-T) identities should be 
imposed on the renormalised theory, and not on the bare quantities (we 
do not consider anomalies). The corresponding renormalised quantities are 
denoted by a sub- or superindex R. We then draw the consequences for the 
specific MOM scheme. 

Z 2 denotes as usual the fermion held or propagator renormalisation ac¬ 
cording to : 


q = s/Z~ 2 qR 

S(p ) = Z 2 S R (p) 


( 8 ) 


Let us recall that the corresponding renormalised vertex functions are defined 
through: 


r(p) = zpZo 1 r E ( P ), 


(9) 


where the necessary subindices are implicit for each type of vertex; Zq is the 
renormalisation of the composite operator, namely a current or density oper¬ 
ator : O = jv 1 jAt P 5 ; the Z 2 factor is to take into account the amputation Q. 

The lattice calculations, being done at a finite cut-off, generate, as other 
regnlarisation schemes, finite 0 (g 2 ) effects, due to additional divergencies 
multiplying the a terms (which have higher dimension), which vanish slowly 
with the inverse cutoff or lattice unit a, and are included in the factors 
Z 2 , Zy, Za- There are also terms with powers of a which we do not write. 
Zy , Za are independent of the renormalisation scheme up to such terms. 
The fact that we do not include such terms means that our equations should 
hold only sufficiently close to the continuum. 

^^Note that the standard definition of renormalisation constants is to divide the bare 
quantity by the renormalisation constant to obtain the renormalised quantity (except for 
photon or gluon vertex renormalisation factors Z\ which we do not use). In principle, 
renormalisation of composite operators, for instance Zy , should be defined similarly. We 
have followed this convention in our works on gluon fields, for the renormalisation of 
A 2 . But, in the case of quark composite operators, an opposite convention has become 
standard in lattice calculations : (qOq)bare = Zq 1 (qOq)n ; we feel compelled to maintain 
this convention for the sake of comparison with parallel works on the lattice. This explains 
our writing of the renormalised vertex function. 
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3.2 The axial Ward identity 

Let us develop the consequences of the axial W-T identity. We define a bare 
mass p through the equation 


d,(.nT = 2 P P 5 


( 10 ) 


(the notation p is old and unsuggestive of a mass, but it avoids any ambiguity 
in a world where there are so many masses). In renormalised form : 


^hUa) r = 2 m R (P 5 ) R 


( 11 ) 


mR, is the renormalised mass in the considered scheme, which then satihes 
the relation : 


- Zp 1 Z A p. 


( 12 ) 


To exploit fully the identities, one has to return first to the general case 
q = p' — p ^ 0. Since they reflect the symmetries of the physical theory, 
the naive Ward identities should a priori hold for the renormalised Green 
functions (except for anomalies) and at infinite cutoff, which means : 


9/x (r A )^(p,q) = -*(^ R 1 (p + g )75 + l 5 S R 1 (p)) + i2m R (T 5 ) R (p,q) (13) 


Returning then to bare quantities which are the ones actually measured on 
the lattice, one gets, multiplying both sides by Z^ 1 and using tur = Z^Z A p 


Z A q fl (r A )»(p,q) = -i(S 1 (p + g)75 + 75-S' 1 (p)) + i 2Z A pT 5 (p, q), (14) 


which depends only on one renormalisation constant Z A , and bare, renor¬ 
malisation scheme independent, quantities. 

Since eqn. (El) has been established without any specification of the 
renormalisation scheme, it shows that Z A also is independent of the renor¬ 
malisation scheme. Therefore, one should expect Z^ IOM to be equal to 
Z A . But this is not the case, as the same identity eqn. (ITT)) shows, see the 
demonstration below. 

3.3 MOM scheme 

Let us now introduce the MOM scheme. It must be first defined for the 
propagator, through conditions at some normalisation momentum p 2 = p 2 , 
originally due to Georgi and Politzer [9]: 



( 15 ) 
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which means 


Zf OM = Z^ix 2 )- 1 (16) 

according to eqn. (JTJ). Also, it means that the renormalised mass is then 

m M° M = (i7) 

i.e. it is the mass function at p 2 = p 2 . 

As to quark current vertices, it is then commonly accepted that they can 
be renormalised analogously by setting (gy\ 5 )r(p 2 — p 2 ) = 1 HI This leads 
to the well-known ”renormalisation constants”, for example : 

= WVs&tP 2 = f) (18) 

The quotation marks are provocative and aim at signalling that they may 
not be the true renormalisation constants. Indeed, it must be stressed that 
having fixed the renormalisation of the propagator, the renormalisation 
conditions of vertices cannot be imposed freely : they must be con¬ 
strained by the renormalised WT identities (similarly to the Slavnov identi¬ 
ties for QCD vertices); one is then not allowed to set (gy\ 5 )r 0 2 = h 2 ) = 1 
freely. 

4 An equation for "Z\ l<)}1 " /Z,\ by derivation 
of the WT identity near <f = 0 

Let us first recall the very simple argument concerning the vector case. In 
the vector case, the Ward identity is very simple : 

Zy <f(IV) M (p, q) = -i (A _1 ((p + q) - S _ 1 (p)). (19) 

Then, as is well known, by derivation with respect to one gets among 
other relations: 


Zy — Zf(p 2 )/gy\p 2 ) ( 20 ) 

where gy\p 2 ) is the coefficient of the 7 ^ term in the Lorentz decomposition 
of the vector vertex. Zy has thus be determined independently of any choice 
of renormalisation scheme : it is indeed independent of the scheme, being 
expressed in terms of bare quantities. 

2 We set standard conditions on one invariant. We are aware that others may be set by 
combining several invariants in a trace. They lead to complications in the discussion. 
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It must be noticed, of course, that the r.li.s. of eqn. (J2Uj) is nothing 
else than the MOM renormalisation constant at p 2 = p 2 dehned by the 
renormalisation condition (gv\g 2 ))ii = 1 , i.e. : 

= “Zf OM (/i 2 )” ( 21 ) 

Therefore, “Z^ OM (/i 2 )” is indeed the expected renormalisation of the vector 
current and we may abandon the quotation marks 0 . 

Therefore also, up to now, everything is well with W-T identities in the 
MOM non perturbative scheme. 

Now, for the axial case, comes the inconsistency . As in the vector current 
case, the axial W-T identity will give a constraint on the axial vertex at q = 0 
by taking the derivative of eqn. [TT] with respect to q at q = 0. We get : 

Za{T a )» = 75 + 2 i Z A p—V b (p, q) (22) 

It must be stressed that not only this relation (\22\i is more complex than 
in the vector case , due to the T 5 (p, q) contribution, but also the latter 
does not vanish in general even in the chiral limit, because of the pion pole 
in T 5 (p,q) a This is one more manifestation of the spontaneous breaking of 
chiral symmetry. Also, of course, it must be recalled that, on the lattice, the 
Ward identity is not exact, but holds only up to artefacts, because we work 
at finite cutoff, and the deviation could be found very large in some cases. 

Let us comment more on the chiral limit. The demonstration has been 
done away from the chiral limit, and the relation is at q^ = 0. Now, the chiral 
limit of the r.h.s. of eqn. (EZZj) is regular since the coupling of a pseudoscalar 
pion to the axial current has a factor q^ : 

(ttICuVIO) oc (23) 

and vanishes at q M = 0. Therefore, the limit of the r.h.s. must be also regular, 
although non zero. This will be shown explicitly below. 

From the equation (l 22 ]h one deduces that Z A ^ where Ll z^ IOM ” 

is dehned, in parallel with Zy IOM , as Zp(p 2 = p 2 )/g^\p 2 = / 1 2 ). This is due 
to the derivative of the pseudoscalar term. In fact Z^(p 2 = p 2 ) /g\\ (p 2 = p 2 ) 
is not even independent of / 1 2 ; one can hope only that it reaches Z A at large 
p ; then, it would be perhaps better to discard this MOM definition, since 
the word is misleading. 

3 It must be observed that we stick strictly to definition of MOM condition through one 
invariant i = 1. We do not consider sums over several invariants as done sometimes 

4 A somewhat different expression was given in the previous paper, due to a confusion 
with an older definition of Z^ IOM through traces. 





Let us show this statement in more detail, by tracking the contributions 
having the 7^75 structure of the undressed vertex, to match the invariant 
g A \p 2 — p 2 ) hi the axial vertex. From now on, we keep to the general p 
instead of setting p 2 = /1 2 , which is not useful. In r 5 (p, q), written in full in 
equation ([ 5 ]) the relevant term is obviously 

r 5 (p, q) = + ••• (24) 

The derivative at q = 0 gives a contribution 2 iZ A p x —ig^ {jp 2 )'y^ to eqn. 
dHH H . Therefore : 

Z A g { l\p 2 ) = Z^p 2 ) + 2 Z A pgf\p 2 ) (25) 

or 

Z A /Zf OM (p 2 ) = 1 + 2 Z A p gf\p 2 )/Z^p 2 ) (26) 

The derivation on the propagator coefficients Z^(p 2 ) and Z^(p 2 )M(p 2 ), as 
well as on the other invariants in T 5 is seen to give contributions to the other 
invariants in the expansion of the axial vertex, eqn. (pi]). 

One can express the result (1261) in a more striking form, returning to the 
renormalised axial vertex, and to g 2 as the MOM renormalisation scale for 
the propagator. Since Z 2 (g 2 ) = 1 /Z^(g 2 ), Z A /Z^ WM (g 2 ) is nothing else 
than Z A Z 2 g A \g 2 ) = (g A > )n ( p 2 = g 2 )- Then the relation is nothing but : 

(g { A )R(p 2 = P 2 ) = 1 + 2 Z A Z 2 (g 2 )p gf\g 2 ) (27) 

which exhibits clearly the statement that (g[l > )R(p 2 = g 2 ) cannot be chosen 
arbitrarily once the propagator has been renormalised : its value is completely 
determined, and in particular it cannot be set to 1 . 

We see no reason why g^\p 2 ) should vanish. Rather, it is clear that it 
contains a pion pole contribution, since 7 M g M 75 is a known structure in the 
Bethe-Salpeter vertex function of the pion, see for instance the appendix of 
Nambu and Jona-Lasinio, Phys.Rev. 124 (1961) 246-254 [7]. It is reason¬ 
able to suppose that the effect vanishes by powers at large p 2 , since this is 
the general behaviour of Green functions. One must note however that the 
transition to 0 as seen from Fig. 5 of j' 6 | is rather abrupt. 

5 One must be aware that since there are other possible definitions of Z ^ OM , involving 
traces, one would obtain for them a similar equation, but with a contribution of different 
invariants of the pseudoscalar vertex, for instance the derivative of ; (p, q) for a trace on 

7a»75- 
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The chiral limit of the deviation is finite but non zero as announced, since 
the pole l/ml of gf\p 2 ) is multiplied by the factor m q %. 

We have not heard of such an inconsistency of the MOM definition of 
the axial vertex renormalisation in perturbative QCD. At least the deviation 
should vanish in the chiral limit in perturbative QCD, since there would no 
longer be a pion pole to compensate the m q factor. On the other hand, it 
remains to be known by explicit calculation what happens at m q 7 ^ 0 . 

Finally, let us note that if one were using definitions of Z\ IOAI with traces, 
as has been usual for some time, one would have to include other invariants 
of r 5 in the expression of Za/Z^ om . 

4.1 Connected observation of momentum dependence 
of Za° M /Z y OM in lattice calculations 

The effect may be exhibited most clearly with chiral invariant actions at 
m q = 0 where Za/Zv = 1 is expected to hold exactly. And indeed, it 
seems to have been seen in certain lattice simulations of Z^ OM jZy OM with 
chiral symmetry preserving actions, showing near the chiral limit a decrease 
from 1 with decreasing q 2 (Dawson, with domain wall fermions [ 8 J ; our 
paper Pliys.Rev.D74:034505,2006 [ 6 J on the quark propagator with Ginsparg- 
Wilson action, especially fig. 5), while it reaches 1 at large momentum. But 
in neither of these works, was it possible to separate cleanly the effect from 
artefacts, and the works should be redone. In connection, a lattice calculation 
of the new invariants in the pseudocalar vertex should be made to ascertain 
the estimate obtained in eqn. (J2HD- 


5 Relation with the discussion of Martinelli 
et al. 

In fact, there is in paper [lj a discussion also concluding to a difference 
gMOM Za- Apparently it has not led to further discussion. It has a 
connection with ours. Their argument also rests on the derivative of the 
W-T identity at q — 0, and they also conclude that the problem should 
disappear at large momenta. However, our finite result does not appear in 
their calculation. 

6 There has been statements that Z^ OM /Zy OM = 1 in the chiral limit, but they rest on 
the assumption that the vacuum is chiral symmetric, so they are valid only asymptotically 
in p 2 


10 



Why ? This requires a rather long explanation. The first difference is 
that they work strictly at m q = 0 , without taking a limit from m q 7 ^ 0 . 
At first sight, this approach is quite opposite and could be expected to be 
incompatible with ours. Moreover, one could feel dangerous to work at m q = 
0 and indeed we have preferred to start from m/0. Nevertheless, we show 
finally that one could obtain the same result as in our approach for m q —y 0 . 
But this requires using the m q = 0 approach differently. 

1) Let us first compare the starting point of their discussion. 

Their argument is based on their eqn. (12) where the W-T identity has 
no pseudoscalar term in the r.h.s. On the contrary, our pseudoscalar 
term in the r.h.s., is non zero even in the chiral limit (in the form 
0/0). This r.h.s. pseudoscalar term is crucial in our argument, and in fact 
the effect is present at any mass m q , not only in the chiral limit. 

Moreover, in their treatment at m q = 0, they find by derivation of the 
W-T identity their eqn. (13), with some non explicited contribution to the 
derivative of g^T^^G 9 , g) at g ^ 1 = 0 due to the q 2 = 0 Goldstone pole in 

(r A) M (p,g)- 

We do not find such a contribution in our treatment with m q 7 ^ 0, because 
(Ta) m G»,<?) is not singular at all at q — 0 , the pole being shifted by m 2 s . 
We obtain rather a g M /(g 2 — m 2 J, therefore a contribution g 2 /(g 2 — m 2 s ) to 
9 M (ru)^(p, q), whose derivative is 0 : 

7^g 2 |gM= 0 = 0 (28) 

(We work here in Minkowski space-time for easiness) 

2) In our opinion, a m q = 0 approach, although perhaps virtually dan¬ 

gerous, is nevertheless possible with various precautions, as we show belowH. 
But there is another problem in [T|. One decomposes a)h into a sum 

of two terms (l.h.s. of their eqn. (13): 

)<• + < 29 ) 

This decomposition has the drawback that both terms of the sum are 
singular at q 2 = 0 , while their sum was regular : the singularities are gM/g 2 
in their first term, —q u l/q 2 in the second one. The finite difference is then 
not made explicit. 

' Let us recall that the mechanism of Nambu-Jona-Lasinio was partly illustrated within 
this strict chiral symmetric situation, at m q = 0 [7) 
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3) Let us now show that one can retrieve our result, at least in the chi¬ 
ral limit, through a calculation strictly done at m q = 0 but avoiding this 
decomposition. 

Let us consider the Goldstone contribution. There is no pseudoscalar 
term in the r.h.s. of W-T identity in the m q = 0 method, and the pole is 
now only present in (r^) M , and it has the form : 

q) (30) 

q 

where T n (p, q) is the full pseudoscalar vertex of the pion i.e. a 4 x 4 with a 
structure parallel to the one of Ts.. 

ft gives a regular contribution to q^{Y a)^ '■ 


q'^f^Trip, q) = fnTnip, <?) (31) 

q 2 

The derivative of q fJj (J'A) fJ , at = 0 then has a term not present at m q ^ 0 : 

fir-^ r n(p, q )\ q v=o (32) 

However, this is exactly equivalent to what we get in our m q ^ 0, —> 0 
method, except that in our case the quasi-Goldstone pole is on the other side 
of the equation Za(X a)h = • ••, see the r.h.s. of our eqn. (l22]h second term, in 
T 5 . The sign should be naturally opposite, since the contribution is on the 
other side of the W-T identity in the m q —> 0 method ; its contribution to 
the pseudoscalar term is: 


( 71 - 1 - 075 - 010 ) 

q 2 — ml 


r n(p,q)\ 


and the contribution to eqn. (l22l is then) : 


(33) 


lim, 


m \—>0 fir 


d tv? d 

rs u ~2 q) |qrM =0 = fir r\ „ rVCP, q) 1^=0 

oqv q 2 — ml oq» 


(34) 


using f n ml = 2Z j 4p(7r|-07 5 '0|O). It is the same as in eqn. (j32]L except that 
the sign is ducly opposite. 

Now, the effect is present at any quark mass but going to the chiral limit 
is much easier by our main method. 
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